The Fate of Bose-Einstein Condensate in the Presence of Spin-Orbit Coupling 
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Intensive theoretical studies have recently predicted that a Bose-Einstein condense exhibits a 
variety of novel properties if spin-orbit coupling is present. However, an unambiguous fact has 
also been pointed out that Rashba coupling destroys a condensate of non-interacting bosons even 
in high dimensions. Therefore, a conceptually important question arises, whether a condensate 
exists in the presence of interaction and a general type of spin-orbit coupling? Here we show that 
interaction qualitatively changes the ground state of bosons under Rashba spin-orbit coupling. Any 
infinitesimal repulsion forces bosons either to condense at one or two momentum states, or to form 
a Schrodinger-Cat state that is a superposition of infinite numbers of fragmented condensates. The 
Schrodinger-Cat is unstable against the anisotropy of spin-orbit coupling in systems with large 
numbers of particles, leading to the revival of a condensate in current experiments. 



Spin-orbit coupling(SOC) is the underlying mecha- 
nism for many fundamental quantum phenomena, rang- 
ing from the atomic fine structure to the newly discovered 
novel properties of topological insulators[TJ [2]. The re- 
cent realization of synthetic SOC for neutral alkali atoms 
in laboratories provides physicists a new platform to 
study the effects of SOC in many-body systems, in which 
a wide range of parameters can be well controlled in 
experiments [3H2]. Particularly, spin-orbit coupled bosons 
offer physicists a unique opportunity to explore how SOC 
may manifest itself at the macroscopic level. As it is 
known as a textbook result that bosons naturally form a 
condensate at the ground state in three and two dimen- 
sions, intensive theoretical efforts have predicted a num- 
ber of macroscopic quantum phenomena exhibited by a 
Bose-Einstein condensate in the presence of SOC[7HT6]. 

On the other hand, a significant effect of SOC on non- 
interacting bosons has also been realized recently. It was 
pointed out that some types of SOC may completely de- 
stroy a non-interacting condensate at zero or any finite 
temperatures even in high dimensions [T7HT9] . To see this 
effect, one can start from the single-particle Hamiltonian 
for spin-orbit coupled bosons, which can be written as 

* = £ / d D r& a (r) [~^Pj *a(r) + H SO C, (1) 

where D is the dimension, a =t?>lo ^3r( r )(^cr(r)) is the 
creation (annihilation) operator at r and M is the mass. 
The part of SOC in the Hamiltonian is described by 

H SO c = -i J d D r^\(r)(\ x d x -i\ y d y )<i> l (v) + c.c., (2) 

where X XjV is the coupling strength. For D = 3 the low 
energy part of single-particle spectrum can be written as 

e k = ^(\k ± \ 2 +k 2 z )- y/xikl + Xykl, (3) 

where k = (k x ,k y ,k z ) is the momentum and = 
(k x , k y ). For D = 2, one can simply set k z = in Eq.([3|. 



Among all the configurations of SOC, Rashba coupling 
corresponding to X x = X y = X is of particular interest. 
In both D = 2 and D = 3, the kinetic energy mini- 
mum under Rashba coupling becomes a circle in the x-y 
plane with radius |kj_| = fco = MX/h 2 , which means an 
infinitely degeneracy of the single-particle ground state. 
Correspondingly, the low-energy Density of States be- 
comes that in (D — 1) dimension in the absence of SOC. 
Therefore, a non-interacting condensate is completely de- 
stroyed at zero temperature for D = 2 [19] and at any 
finite temperature for D = 3 [T7( [T8]. 

Whereas the disappearance of a non-interacting con- 
densate with Rashba coupling is rather clear, one may 
naturally ask whether interaction will change the fate of 
a condensate. To answer this question, we first note that 
the infinite degeneracy of single-particle ground states 
makes interaction effects highly non-pert urbative. To 
determine the many-body ground state amounts to se- 
lecting the state with the lowest interaction energy from 
an infinitely degenerate subspace, which is known to be 
a challenging problem. In the literature, a mean field ap- 
proach has been adopted to compare energies of a special 
class of states that can be described by condensate wave 
functions [8- 14 . However, with the observation that SOC 
may completely destroy a condensate, one naturally con- 
cerns the validity of mean field approach. To explore the 
exact many-body ground state, interaction effect needs 
to be investigated more systematically. 

We expand the interaction between two-component 
bosons 

U = n J d D r( gi n 2 t (r) + 2 nf (r) + # 12 n t (r)%(r)) (4) 

in the basis of single-particle eigenstates, where h a (r) = 
^J.(r)^ CT (r) is the density operator, Q is the volume of 
the system, and ^=1,2 > 0, gi2 > are the intra- and 
inter-spin repulsion. In this basis, the effective interac- 
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FIG. 1: Hartree-Fock energy E^l of N particles. Red circle represents the kinetic energy minimum, and blue ellipses represent 
a large number of bosons occupying the same momentum state. Blue and white dot represent a single boson and a hole in the 
condensate respectively. Single-particle and fragmented condensate minimize E^/p for 7 < 1 and 7 > 1. They are separated 
from high energies states by an energy gap Eg = g + N5g/2 or Eg = g(N/2 — l), where g = (g± +#2 +#12) and Sg = gi +#2 — #12- 
Any two fragmented condensates with even particle numbers N — 2n are coupled through a sequence of tunneling induced by 
n steps of scattering, as shown by the green arrows. 
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(L 9 ) is the creation (annihilation) operator in the 
subspace C composed by all single-particle ground states, 
and (j) is the polar angle on the circle |kj_| = fen [19^ [20] . 

Using Eq.([5|, a number of significant interaction effects 
can be revealed. First, any infinitesimal repulsion forces 
bosons to form a single-particle condensate or a frag- 
mented condensate, distinct from non-interacting sys- 
tems where a condensate is absent. Second, interaction 
inevitably builds up a superposition of fragmented con- 
densates with even particle number, and the exact ground 
state becomes a Schro dinger- Cat beyond the prediction 
of mean field theory. Finally, Schrodinger-Cat states col- 
lapse in systems with large numbers of particle if a finite 
anisotropy exists in SOC, leading to the revival of a con- 
densate in current experiments. 

We start from the angular dependent Ul in Eq.^, 
which drastically changes the Hartree-Fock energy of in- 
teracting systems. To see this fact, a two-body prob- 
lem is an illuminating example, where Fock states can 



be written as either |0, 6) = L^Z/^ +6 ,|0) for ^ or 

|0,O) = L^LjjO) /\/2. 6 is the relative angle between 
the two bosons on the circle. The Hartree-Fock energy 
E 
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the result for a TV-body problem, can be written as 
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Eq.|7| shows that, |0, 0) and |0,7r) minimize E^ F for 
7 < 1 and 7 > 1 respectively, where 7 = #12 /(#i + #2)- 
States with ^ 0,7r always have higher energies. For 
7 > 1, two bosons occupy two opposite points on the cir- 
cle, i.e., fragmentation occurs. Without SOC, it is well- 
known that fragmentation is not possible in homogenous 
systems, because of the cost of Fock energy when bosons 
occupy different momentum states [21]. The presence of 
Rashba coupling suppresses the term in Fock energy that 
is proportional to g\2 by a factor of cos 0. Therefore, frag- 
mentation may occur at large enough #12. 

The above discussions can be directly generalized to a 
TV-body system. For systems with even particle number, 
a single-particle condensate 

\c^u = mo)/m, 



and a fragmented condensate 



\F^U = L^mZ 2 \0)/(N/2)\ 



(8) 



(9) 
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minimize E^ F for 7 > 1 and 7 < 1 respectively, where 
<p is arbitrary, reflecting the rotation symmetry in the 
momentum space for Rashba coupling. This result can 
be understood from that any pair of bosons prefer to 
occupy the same or the opposite points on the circle to 
save interaction energy. For odd particle number TV = 
2n + 1, where n is an integer, the same conclusion holds, 
with the expression for the fragmented condensate simply 
replaced by 



Lt"Lt" + V|0)/v^ 



+ 1)!. (10) 



Any other Fock state in C costs additional Hartree- 
Fock energy proportional to the total particle number, as 
shown in Fig.l. Therefore, all bosons condense at one or 
two momentum states, distinct from non-interacting sys- 
tems where bosons can distribute on the circle |kj_| = ko 
arbitrarily. 

While the significance of interaction effects has already 
been reflected in Hartree-Fock energy, it is worthwhile to 
point out that E^ F only takes into account the diagonal 
term of interaction matrix elements. For single-particle 
condensates and fragmented condensate with odd par- 
ticle number, such a consideration is sufficient. States 
|CM)0 or |F[ 2n+1 ]) with different values of carry dif- 
ferent total momentum. Therefore, the minimums of 
Hartree-Fock energy cannot be mixed with each other 
by interaction due to constraint of momentum conser- 
vation. However, such a constraint is absent for frag- 
mented condensates with even particle numbers, as all 
states \F^ 2n ^) ( f ) have zero total momentum, and any two 
of them are inevitably coupled by a series of scattering, 
as demonstrated in Figure 2. This type of tunneling is 
known to be crucial in the presence of degenerate ground 
states, and produces macroscopic or mesoscopic quan- 
tum coherence in various systems [22- 24 . In our case, it 
fundamentally changes the ground state structure when 
fragmentation occurs with even particle number. 

We derive an effective Hamiltonian in the subspace 
composed by all fragmented condensates (See Appendix), 



H eff = I #4" nJ olO0+ / ##%,^6I,6*(i-(W)> 
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(11) 

where = L^L^^/n! is the creation operator for a 
fragmented condensate, and 
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The diagonal term S F n ^ = E^p + where E^p is the 
Hartree-Fock energy of a fragmented condensate, and £ ' 
is the diagonal part of the correction to the energy due to 
the coupling between \F^ 2rx ^)^ and high energy states(See 
supplementary materials). As 6' is independent on </>, it 
does not affect the structure of ground state wave func- 
tion. 



While H [ ?f) is analytically solvable, the underlying 
physics can be revealed by considering the simplest case 
with #1=0, where the system can be mapped to a one- 
dimension ring with an "infinite-range" constant tunnel- 
ing t = (-l) n_1 #2 n2 /(#2 + #i 2 ) n-1 [ 2 5 • The eigenstates 
can be rewritten as 



\SC^} = e i2 ^O+|0), 



(13) 







where v is an integer, with the corresponding eigenenergy 

42 „] = e p»] + { _ ir -i _^£_ { ^ uo _ 1} _ (14) 

As the reduced single-particle density matrix (L^L^,) of 

is zero, and the reduced 2n-particle density ma- 
trix (L^ n L^ 7r L^ / L^ /+7r ) has a unique macroscopic eigen- 
value (n!/7r) 2 , \SC [ u n] ) is a Schrodinger-Cat state, which 
has distinct properties from ordinary single-particle or 
fragmented condensates [24]. For general cases where 
both gi and #2 are finite, the structure of Schrodinger- 
Cat is similar but more elaborated(See Supplementary 
Materials). 

It is known that Schrodinger-Cat state is fragile in 
systems with large numbers of particles. In our case, 
this very nature of Schrodinger-Cat is directly related 
to the fate of a condensate in the presence of SOC. 



Eq.(14) shows that the energy difference between the 
Schrodinger-Cat state and the fragmented condensate 
can be characterized by 77 = 1 — E^ n ^ /£p U \ For a small 
total particle number N = 2n, 77 takes a large value. 
For example, one estimates that 77 ~ 0.61,0.11,0.1 for 
N = 4, 6, 8, and #2/gi2 = 0.9, which means a large energy 
gained by forming a Schrodinger-Cat. While this fact 
strikingly changes the properties of a few-body system 
(See Supplementary material), one notes that 77 decreases 
exponentially with increasing particle number N = 2n, 

/ \-n+l 

i.e., r] ~ 1 1 -f ^ \ . This is a typical feature of a 
high order process that requries multiple steps of two- 
body scattering to couple two degenerate states. When 
N — >• 00, 77 vanishes, which means a fragmented con- 
densate becomes degenerate with the exact many-body 
ground state in the thermodynamic limit. In current ex- 
periments, N ~ 10 5 — 10 6 , one can easily see that 77 is a 
negligible number. 

The tiny 77 for a large value of N indicates that a big 
Schrodinger-Cat is unstable against external perturba- 
tion. For our discussions, the anisotropy of SOC, which 
always exists in practice, is naturally such a perturba- 
tion. Assuming > A y , Eq.(|3| shows that kinetic 
energy minimums are located at two separate points 
k = ±k x o = ±MX x /h 2 . We define the circle |k^| = 
k x o as £ ', on which the expression for interaction en- 
ergy is identical to that in Eq.(ll), with <j) replaced by 
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(p = 8iYg{\ x k x , Xyk y }^9\. The kinetic energy in CJ is now 
a function of <b, 



K{4>) = m & 



k x0 X x Jl-a(2-a)sm <j>, (15) 



where a = (X x — X y )/X x characterizes the anisotropy of 
SOC, and leads to an offset of kinetic energy A(a) = 
(K(n/2) - K(0))N = aNMX 2 x /h 2 in 

For single-particle condensates and fragmented con- 
densates with odd particle numbers, any infinitesimal 
a picks up | (7^)^=0 or |F^)^ = o as the ground state, 
which favors the kinetic energy with no cost of inter- 
action energy. For Schrodinger-Cat states, we define a 
characteristic anisotropy, 
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(16) 



by setting the kinetic energy offset equal to the 
strength of interaction induced off-diagonal coupling, i.e., 
A(a*) = t. For small anisotropy a <C a*, interac- 
tion energy is dominate and Schrodinger-Cat forms. For 
a ^> a*, to form a superposition of fragmented con- 
densates costs too much kinetic energy, and a single 
fragmented condensate |F 1^)0=0 becomes the ground 
state. A more quantitative analysis is given in Supple- 
mentary Materials. Using Eq.([l6|), one estimates that 
a* = 9%, 4%, 2% for N = 4,6,8, g 2 /gi2 = 0.9 and 
g 2 nh 2 /(2MXl) = 0.2. This means a small Schrodinger- 
Cat is stable in a finite range of anisotropy. With increas- 
ing N, a* decreases exponentially. One can verify that 
for N ~ 10 5 — 10 6 , the typical particle number in current 
experiments, a* is essentially zero. Therefore, one con- 
clude that a Schrodinger-Cat collapses to a fragmented 
condensate. An additional external potential that breaks 
the translation invariance shall further change the frag- 
mented condensate to a stripe condensate(See Supple- 
mental Materials). 

Whereas we have answered fundamental questions on 
whether, when and why a condensate exists in the pres- 
ence of SOC, our work also demonstrates that the cou- 
pling between Hartree-Fock energy minimums becomes 
particularly important in a spin-orbit coupled few-body 
system. Evidence for this type of coupling has been found 
in a recent numerical calculation [215]. Provided that sig- 
nificant progresses have been made on manipulating a 
few atoms in current experiments [27j [28], we hope that 
our work will stimulate more studies on SOC induced 
novel ground states of mesoscopic cold atomic systems. 
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Physics, The Chinese University of Hong Kong. X.C. 
acknowledges support from the Initiative Scientific Re- 
search Program of Tsinghua University and NSFC under 
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Appendix 



We define the excited state with n — m pairs of bosons 
at the angle <p and m pairs of bosons at <j>' as \m) ^ . 
Within this definition, \F^}^ = |0)^/, |F [2n] )/ = 
|n)<^0/. Interaction scatters a pair of bosons out from 
|O)0, 0' to the momentum states at k 7 = (/co,</>') and 
— k 7 = (fco,^' + 7r). The Schrodinger equation satisfied 
by lO)^^/ can be written as 



E)\0)< 



#Uo|lW =0, (17) 



where U mn = </>,</>' (ra|Wz,|n) , and E is the eigen en- 
ergy. Considering the shorted path that connects |O)<^0/ 
and \n) ( f )y( f ) / in Hilbert space as shown in Fig.(l), the 
Schrodinger equation satisfied by {m)^' can be approx- 
imated by 



(Em - E)\m)^^ + ^2 Um+i,m\m + 00,0' = °> ( 18 ) 
l=-l 

where E m = {m\U\m) is the Hartree-Fock energy 
of |m)0 } 0/. Eliminating \m) ( f ) ^ from the above equations, 



the effective Hamiltonian in Eq.(ll) that takes the stan 



dard form of a high order perturbation theory is obtained. 
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SUPPLEMENTARY MATERIAL 

In this supplementary material, we present the results 
on the diagonal term of Eq.(4) in the main text, the gen- 
eral solution for the effective Hamiltonian H e ff, small TV 
systems, effects of anisotropy of spin-orbit coupling, and 
the evolution from a fragmented condensate to a stripe 
condensate. 

Diagonal term of Eq.(4) 

In Eq.(4) of the main text, the Hartree-Fock energy 
can be written as 

e hf = \(6gi+6g2 + 2g 12 )n 2 -^( gi +g 2 +g 12 )n, (19) 

For £ to the second order approximation for the diago- 
nal term, we consider the lowest excited states with only 
one pair of atoms scattered out from the fragmented con- 
densate and obtain 

£' = C d4>'\(Ex\U L \F^)\ 2 /E G , (20) 

JO 

where \Ex) = L^ 1 h\, h\, + J0) / \n - 1)! and 
E G = (Ex\U L \Ex) - ={ gi + g 2 + <7i 2 )(n - 1). To 
be explicit, £' = tt^ 3 ' +g " 



9i+92+9i2 ' 



Solution of H e ff 

For general cases of gi ^ 0, gi ^ 0, V^p can be written 



n-l„2 n 



Vd>.d>' — 



52+512) 



-Y, C n9i92-^ 2i ^'-^. (21) 
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FIG. 2: Many-body ground states under Rashba coupling, 
single-particle condensate occupies the region with small 7. 
For large 7, the ground state is a fragmented condensate at 
odd particle number, TV — 41 + 1 or TV = 41 + 3, where / 
is an integer. For even particle numbers, the ground state 
is a Schrodinger-Cat state with infinite degeneracy for (I) 
with TV = 41 + 2 and no degeneracy for (II) with N — 41. 
The dashed red curve represents the critical value 7 C for the 
transition between single-particle condensate and other three 
states. 



It is straightforward to verify that \SC\? n ^) are the eigen 
states with the eigen energies 



£ [2n] ^ (-1)"-^ 



E [2r t 



F - (./^^^ (^n^^^ - (91 + 92T) 

< v < n 

p [2n] _ / -.\ n -l {gi+92) n n 2 
°F V L ) (g 2 +gi2) n - 1 

v < 0, v > n 



K = 



(22) 

where C v n = v ^_ v)v For TV = 2n = 41 + 2, (-l) 71 ' 1 > 0, 
all states with v > n are degenerate and have lower 
energy than those with v < n. This means the ground 
states are infinitely degenerate. For N = 2n = 41, 
( — l) n_1 < 0, the ground state is unique, which is 
determined by the ratio of g\/g2 so that CX l g\g^~ v 
is maximized. For example, for g\jg 2 = 0,1, 00, the 
winding number for the ground state = 0,n/2,n 
respectively. 

Small TV systems 

The dependence of the ground state on the total par- 
ticle number and interaction is summarized in Fig. (2). 
7 C (TV) is the boundary between different phases, which is 
determined by comparing the energy of the single-particle 
condensate \C^) with other states. For instance, by 
solving E% n] = Ej? n \ one obtains the value of 7 C (TV) 
that separates \C^) and the Schrodinger-Cat \SC^). 

At small particle numbers, the derivation of 7 C (TV) 
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from the mean field value 7J? = 1 is significant. For TV = 
2n = 2, this two-body problem can be solved exactly. Un- 



der this situation, £' is absent, and £ F 



[2] 



E HF = 91+92- 

By setting n = 1 in Eq. p2j) , it can be verified easily that 



for any v 7^ 0,1, the eigen energy is exactly zero. Com- 
paring with the energy of any state with a finite total mo- 
mentum, we see that the ground state is always a small 
Schro dinger- Cat state, regardless of the ratio between g\2 
and #1 +# 2 . 

In spite of the repulsive interaction, we see that the 
ground state energy of two bosons for Rashba coupling 
is exactly zero. The underlying physics is that the scat- 
tering of a pair of bosons with zero momentum on the 
circle cancels the positive Hatree-Fock energy. It is in- 
teresting to note that, this result is consistent with the 
conclusion of a completely different approach based on 
renormalization of two-body interaction, which discov- 
ered that the effective interaction between two particles 
on the opposite points of the circle is zero [20]. 

It is also useful to study how 7 c (iV) approaches the 
mean field value at large N. Clearly, the contribution 
from £ ' to the ground state energy is much larger than 
that from off-diagonal term in Eq.(4) of the main text. 
If one ignores £' \ E^ = £ F n ^ leads to the mean- field 
result 7^ = 1. Taking into account £' \ the value of 7 C is 
modified. In large N limit, we find 



7 C (2n) « 7; 



.0 71 91+92 



n (9i+92) 



2 ' 



(23) 



Similar scaling of 7 C (A/") — 7^ ~ 1/N holds for odd 
particle numbers. 

Effect of anisotropy of SOC 

The kinetic energy offset in C as shown in Eq. (8) of 
the main text leads to an additional term in the Hamil- 
tonian for anisotropic SOC, 



/* 
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d<j>K{<p)L\L^. 



(24) 



The matrix elements of W between | SCl? ni ) can be cal- 
culated straightforwardly, 



n r 

7T Jo 



d<pe 



2i{v-v' 



(25) 



Since K(4> + tt) = K(cj)), W v y is finite only when 
v — v' is an integer. For Rashba coupling, K((f>) is 
a constant and therefore W u y = 0. With increasing 
5K = K(tt/2) — K(Q) or equivalently the anisotropy of 
SOC a, the amplitude of W u y increases. 



Tthe general expression of the ground state wave func- 
tion can be written as 



\SC™) a = J #<^6J |0>. (26) 

For isotropic SOC and N = 2n = 4/, C = e 2 ™ * / y/n , 
and I 1 2 = 1/tt is a constant, where v$ < n as discussed 
before. In the presence of a small anisotropy, VV is fi- 
nite, one may take into account only the mixing between 
and l^cj^^), which leads to a small variance 
of IC^I 2 . With increasing a, more states with different 
v are mixed into the ground state. As W v y ~ NSK, 
we see for large enough NSK ^> E v>n — E VQ1 all states 
|5Cj> 2n ') are mixed. An equal superposition of them leads 
to C<f) = ^0,o/aA" an d the ground state becomes a frag- 
mented condensate. 

For N = 2n = 41 + 2, a subtle difference emergies due 
to the infinitely degenerate ground states under Rashba 
coupling. In the presence of any infinitesimal VV, a su- 
perposition of these degenerate state \SC^>\) may fulfill 
the goal of minimizing both the kinetic energy and inter- 
action energy. For instance, in the simplest case where 
gi = 0, the ground state becomes 



\ SC [2n] 



CMO), (27) 



where C<j> = (ttS^o — V)j\/i\ 2 — tt. With a finite a, 
\SC^) will be mixed to the ground state. For large 

- E UQ , similar to the case of 



enough NSK > E v 



N 



2n 



4/, all states \SC [ „ n] ) are mixed and a 



fragmented condensate becomes the ground state. 

Fragmented condensate to a stripe condensate 

The fragmened condensate is a Fock state in the mo- 
mentum space, with half particles occupying the state 
k and the other half occupying — k . In homogenous 
systems, these two states can not be mixed due to mo- 
mentum conservation. If an external potential exists and 
translational invariance is broken, an additional term ex- 
ists in the Hamiltonian H' = t'L^L-^ -\-c.c, where t' is 
proportional to the strength of the potential. Therefore, 
it becomes a double- well problem in momentum space. It 
is well known that the tunneling between the two "wells" 
is enhanced with increasing particle number N due to 
the bosonic enhancement effect, i.e., H' ~ t' N . In the 
thermodynamic limit, any infinitesimal t' leads to the for- 
mation of a coherent state ~ (Z^ + L^_^ ) N \0) , so called 
stripe condensate in the literature. 



